function [xep,xop]=circevod(x)
    if any(imag(x) ~= 0)
        error('x非实数序列')
    end
    N=length(x);n=0:(N-1);
    xep=0.5*(x+x(mod(-n,N)+1));
    xop=0.5*(x-x(mod(-n,N)+1));
end